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We investigate the interaction between an infinite cylinder and a free fluid-fluid interface governed 
only by its surface tension. We study the deformation of an initially flat interface when it is deformed 
by the presence of a cylindrical object, tilted at an arbitrary angle, that the interface "totally wets". 
Our simulations predict all significant quantities such as the interface shape, the position of the 
contact line, and the force exerted by the interface on the cylinder. These results are compared with 
an experimental study of the penetration of a soap film by a cylindrical liquid jet. This dynamic 
situation exhibits all the characteristics of a totally wetting interface. We show that whatever the 
inclination, the force is always perpendicular to the plane of the interface, and its amplitude diverges 
as the inclination angle increases. Such results should bring new insights in both fluid and solid 
mechanics, from animal locomotion to surface micro-processing. 
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INTRODUCTION 

Free fluid-fluid interfaces produce fascinating patterns 
in nature. From the simple geometry of drops, bubbles 
and micelles [l| , to the complex but regular arrangement 
of films inside a liquid foam [2| , the concept of interfacial 
energy allows us to rationalize the change of bulk prop- 
erties at the boundaries between two immiscible fluids. 
In the case of liquid-gas interfaces, deformation of the 
interface by a solid object and the subsequent measure- 
ment of the force exterted by one on the other is still a 
widely-used technique to quantify the surface tension (for 
instance the Du Noiiy Ring and Wilhelmy plate methods 
[3|). In many cases the object that pierces the interface 
is cylindrical, and the deformation of a flat liquid-gas in- 
terface due to the contact of a solid rod is a situation of 
interest in both solid and fluid mechanics. 

The propulsion of insects that walk on water is pos- 
sible because of the force between the interface and the 
rigid cylindrical legs [^-lZI- If the latter are flexible, then 
a new class of so-called "elasto-capillary" or "wet-hair" 
problems |8| arise with major applications in nano- or 
microdevices. Thus hairy surfaces, such as carbon nan- 
otubes and biological filaments, are examples of the rich 
morphologies that arise from individual hair deformation 
and collapse [9 4l2 and of the patterns seen in assemblies 
of hairs [13l4l6l | in the presence of a meniscus. Related 
problems are encountered in surface processing, for ex- 
ample in fiber coating and in describing the wetting and 
cleaning properties of textiles 17H19J ; or more specifically 
in probing the wetting properties of nanoscalc droplets 
by deforming them with a cylindrical AFM tip [19|, |20| . 
Finally, the mutual force between a soap film and a cylin- 
drical liquid object, namely a micro-jet, was also found 
to be strong enough to cause the jet to deviate and to 
control its trajectory [21|. 



In general, the rise of a meniscus is fast [221 compared 
to other relevant time-scales (solvent evaporation, leg 
motion). Consequently, the equilibrium shape of the gas- 
liquid interface and the minimum-energy principle are 
usually used to account for observations. 

Here, we use Surface Evolver [23| simulations to study 
the deformation of an initially flat interface when it is 
pierced with an infinite cylinder. We report the effect 
of the inclination angle on the meniscus shape, the con- 
tact line position and on the force generated by the in- 
teraction. In the following we consider the case of a 
"total wetting" condition, meaning that the contact an- 
gle between the interface and the cylinder is zero. The 
results are general as long as the interface is governed 
only by surface tension, meaning that they can be di- 
rectly applied to interfaces between two immiscible liq- 
uids, boundaries between a liquid and a gas, and to most 
2D dimensional lamellae, including soap films, Langmuir 
monolayers, lipid bilayers and biological membranes. The 
main control parameter is the inclination of the cylinder 
with respect to the plane of the unperturbed interface. 
The numerical study is compared with experiments in 
which a cylinder of liquid is brought into contact with a 
soap film; above a certain liquid velocity the results are 
in good agreement with the model, offering further in- 
sight into the dynamical interaction between a cylinder 
and a film. As explained below, the more the cylinder 
is inclined, the greater the deformation of the meniscus 
and the greater the force exerted by the meniscus on the 
cylinder. This force is always perpendicular to the inter- 
face and its amplitude increases rapidly as the inclination 
is increased, diverging as the angle approaches 90°. 

The manuscript is organized as follows: the numerical 
and experimental systems are first presented, then we 
describe and characterize the position of the contact line 
in the simulations, compare with the experiments, and 



a) Simulations 
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sures that the contact angle is slightly greater than zero, 
alleviating numerical problems. 

The surface is meshed with many small triangular 
facets. We tested two levels of refinement, the first with 
about 12,000 triangles, and the second with an extra 2000 
triangles in the region where the rod pierces the sur- 
face; the refinement also controls the number of points 
around the contact line, giving 46 and 81 points respec- 
tively. Both triangulations give equivalent results (data 
not shown), which validates the robustness of the ap- 
proach, and we show the more refined results below. 

One iteration consists of tilting the rod by 3° in the x — 
z plane and finding the new minimum area configuration 
of the surface (Fig. [T^). At each iteration, the location of 
the points along the contact line is recorded which, as we 
shall show below, is sufficient to calculate all quantities 
of interest. 




FIG. 1. Side view showing how the interface deforms, a) 
Simulation, with 9i = 30°, and an imposed contact angle of 
90° (top) or 0° (bottom), b) Experiment, with 6*, = 30°, 
R = 200/im and V = 2m.s-\ 



show how to calculate the force on the cylinder from the 
position of the contact line. We determine the force as 
a function of the inclination from the simulation results, 
and discuss the success of different models in predicting 
this dependence. 



MATERIALS AND METHODS 
Simulations 

We use the Surface Evolver [23[ to calculate the shape 
of a single area-minimizing surface that is pierced by a 
rod at an inclination of Oi to the vertical. All units are 
dimensionless, without loss of generality. We first create 
a planar surface in the x — y plane with surface tension 
1.005 in a fixed square frame of side-length 5.0 and a 
cylindrical rod, with axis in the vertical (z) direction and 
radius R = 0.1, through its centre. To ensure that the 
rod is wetted, we include a virtual liquid film on the lower 
half of the rod with unit tension. The small difference 
in tension between the surface and this virtual film en- 



Experiments 

To counteract the effects of bouyancy (which domi- 
nate capillary forces at these length-scales), we consider 
a double interface formed by a soap film, with surface 
tension 27, rather than the single interface separating a 
liquid from air. Our conclusions remain valid for such an 
interface, but with half the value of the surface tension. 

The force exerted on a static glass rod by a meniscus, 
with a contact angle equal to 90°, is zero. To generate 
a contact angle of 0°, we therefore use a laminar liquid 
jet directed through a soap film of the same composi- 
tion, as described in 2l|, |2J| . The flow in the jet plays a 
subtle role: the liquid flow within the soap film is slow, 
so the film adopts a quasistatic equilibrium shape, yet 
the boundary conditions are governed by the dynamics 
of the jet. As wc shall demonstrate below, friction causes 
the contact line to move downstream, and in the limiting 
case imposes the required contact angle on the soap film. 

The liquid is a soap solution obtained by adding 5% of 
commercial dish-washing liquid (Dreft, Procter & Gam- 
ble) to deionised water. The surface tension of the soap 
solution is 7 = 26.2 ± 0.2 mN.m~^ and its density is 
p = lO'^kg.m^'^. As reported elsewhere [21|, this situ- 
ation is stable and the meniscus forms an equilibrium 
shape which is recorded by a digital camera (Fig. [TJj). 
The incident jet is characterized by its incident angle 0^, 
velocity V and radius R. 

The Weber number, We = pV'^R/"f is used to estimate 
the relative contributions of inertia and capillary forces. 
Within our experimental parameter range, the Reynolds 
number based on the jet characteristics is always signifi- 
cantly larger than unity. 




FIG. 2. Schematic representation of the system. 



RESULTS 



0.15 

0.1 

0.05 



-0.05 

-0.1 




-0.45 
-0.5 

-0.55 
-0.6 

-0.65 
-0.7 





-0.1 0.1 






x' 




0.45 


/-""X 


-0.45 


-0.5 


/ \ 


-0.5 


0.55 
-0.6 
0.65 


M,'.' + oj 


^ -0.55 

-0.6 

-0.65 


-0.7 




-0.7 






0.1 




Shift of the contact line 

The shape and position of the contact hne where the 
interface meets the cyhnder is sketched in Fig. [2] The 
reference frame of the initiahy flat interface is TZ {Oxyz), 
with Ox horizontal, and the reference frame of the cyhn- 
der is TV {Ox'y'z'); they are related by 
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sa\{Oi)x -I- cos(6'i)z. 



(la) 
(lb) 
(Ic) 



All calculations are performed in the reference frame of 
the cylinder, 7?.', using cylindrical coordinates (r', Q' , z'), 
so that any point M has OM = x'ex' + y'cj,' -I- z'cz' — 
r'ev' + O'egi + z'Sz', with 



x' = r'cos{9') 


(2a) 


y' = r'sm{9') 


(2b) 


Cr' = cos(6l')er' + sin(6'')ey/ 


(2c) 


eg' = -sin(6i')4'+cos(6l')ey'. 


(2d) 


y, the contact line (Fig. ^ is given by 




r' ^ R 


(3a) 


z' = Ae'i 


(3b) 



with 0' e [-7r,7r]. 

The position of the contact line in the simulation with 
an incident angle of 9i = 30°, referred to 7?,', is shown in 
Fig. [3] The deviation of the contact line is defined by 



0). Ax'y', Ax'z' and A. 



■y'z' 



Az' 

denote the absolute values of the area enclosed by the 
contact line when projected on the {Ox'y'), {Ox'z') and 
(Oj/'z') planes respectively. By symmetry, Ax'y' ~ 7ri?^ 
and Ax'z' = 0. Ay'z' and Az' are equal to zero when 



FIG. 3. The contact line projected on different planes for a 
simulation with for 6i = 30°, in the reference frame of the 
cylinder TZ' . 



6i ~ and both increase monotonically as 9i increases. 
Fig. Ulshows how Az' /R depends on 9i in the simulation 
(dashed line). 

Side views of the meniscus were used to investigate ex- 
perimentally the effect of R, 9i and V on the meniscus 
shape and perform measurements of Az' . For a given 
set of parameters (i?, 9i), Az' was found to decrease as 
V increases (insert to Fig. |4]) until saturating above a 
limiting velocity. This marks the onset of the dynam- 
ical "total wetting" condition (see Discussion). In this 
regime, the meniscus shape is found to be independent 
of the jet velocity and the following measurements were 
performed at sufficiently high velocity to ensure that Az' 
reached its saturation value. 

Fig. Ulshows Az' jR as a function of 9i, illustrating the 
agreement between simulations and experiments. This 
scaling allows the collapse of all experimental measure- 
ments (several R and 9i values) on the curve found from 
the simulation. 



Force measurement 

Assuming total wetting, it is possible to infer the force 
on the cylinder from the contact line position. 

We denote by M(i?, 9' . z{9')) a point on the contact 
line and dM an infinitesimal variation of position along 
the contact hne (Fig. [2]): dM = Rd9'ee' + dz'cz' = 
d9' [Regi + ^ez' ) . Knowing that the local force is per- 
pendicular to dM and to eV' (total wetting condition). 
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FIG. 4. The quantity Az' /R versus 9i in simulations (dashed 
hne) and in experiments (different radii R — 140, 200 and 270 
fim are shown as circles, triangles and squares respectively). 
Insert: experimental data for Az' versus the incident velocity 
V for a radius of 270 fim and two incident angles: 39° (filled 
circles) and 49° (open circles). The two solid lines indicate the 
values to which the curves saturate. 
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FIG. 5. The different force components Fx and Fz are nor- 
malised by 2'yTTR. Simulations are shown as red stars. Three 
models are tested for F^. Black solid line: vertical force 
model, Fz — 2^ttR/ cos{di). Blue dashed line: total perimeter 
model. Cyan dashed-dotted line: flat interface model. 



the local force dF applied to dM is 

dF ^ -fCr' A dM. 



(4) 



The prefactor 7 holds for single interfaces and should be 
replaced by 27 without loss of generality for double (e.g. 
air-liquid-air) interfaces such as soap films. The total 
force that the interface exerts on the cylinder is 



F 



7 a) eV' A dM, 



(5) 



where C represents the oriented contour around the con- 
tact line. After some manipulation, described in the Ap- 
pendix, an exact expression for the force in the reference 
frame TZ of the interface is obtained: 



F.^ = 2'YTTR cos 



A. 



y'z' 



Fy = 0, 



F. = 277ri? sin(6',) 



27ri?2 
'2^ 



cos 



(6a) 
(6b) 
) ) , (6c) 



where Ay'z' is the absolute value of the projected area of 
the contour in the (Oy' z') plane. Force measurements 
are then obtained numerically using the contour profiles 
calculated in the simulations and the results are shown 
in Fig. [5l For the whole range of 9i tested, F^ is small 
compared to Fz , suggesting that the force acts vertically 
(in the reference frame of the interface) . 

The vertical component of force can be compared to 
three different models (Fig. [5]): 



vertical force model: from the observation above, 
we might assume F^ ^ i). Consequently Eq. [6al 

A I I 

tan(0,) and thus 



gives -^^ 



i^, =27^i?/cos(0,)- 



(7) 



• total perimeter model: the length of the contact 
line, "P, which is measured in the simulations, might 
determine the vertical force: F^ = ^V . 

• flat interface model: F^ could depend on the 
perimeter of the ellipse that results from pro- 
jecting the cross-section of the cylinder on to 
the unperturbed interface {(Oxy) plane): F^ = 
I'P ellipse.- This cUipse has semi-axes of length 
R and R/ cos(6i). An approximation to the 
perimeter of this ellipse (which is valid to better 
than 5% in our parameter range) leads to Fz ~ 

2jnR (l + - (^7^ - 1 
\ TT vcos(6'j) 



DISCUSSION 
Comparison of models 

The deformation of an interface by an object of char- 
acteristic size L should result in a force that scales as 
7L. A rod has two well-separated characteristic lengths 
(radius and axial length) so the force could, in principle, 
range over several orders of magnitude. In the case of an 



infinite rod, one of the characteristic lengths is infinite, 
which explains why the force diverges as the rod becomes 
horizontal. 

It is clear that the vertical force model gives the best 
fit, and the good agreement with the measurements ex- 
tends up to very large inclination angles. The total 
perimeter model is plotted as it represents an upper 
bound for F^, since it assumes that all local forces dF 
(Eq. m are vertical in the reference frame of the in- 
terface. The fiat interface model underestimates the 
force measurements. A priori we might have expected 
the opposite, as for the total perimeter model (which 
also accounts for vertical local forces), but it signifi- 
cantly underestimates the real perimeter of the contact 
line: the presence of the cylinder perturbs the contour 
and there is a significant shift in the vertical direction 
Az = z{9' = it) — z{9' = 0) which increases its perime- 
ter. The vertical force model works well, and it would 
be interesting to see if this is also the case for "partially 
wetted" contact. 

Based on an energy balance |9| a vertical force model 
has also been used in the study of an elastic plate 
that pierces an interface, enabling the authors to infer 
a l/cos(6'i) dependence to explain the different buck- 
ling morphologies of the plate under the effect of cap- 
illary forces. Expression ([7]) could also be helpful in un- 
derstanding the detachment of liquid droplets from fi- 
bres |17Hl9j . in which a drop is pierced by a rod non- 
axisymmetrically. This complicated geometry has two 
contact lines that prevent the drop from being detached 
by gravity, and the rod is locally inclined to the surface of 
the drop. No exact expressions for the interface shape or 
contact line position are available, so that in [17| an ex- 
pression for the force was based on experimental data, 
with no divergence as the inclination angle increased. 
This is in contrast to our 1/ cos{9i) dependence for an in- 
finite planar interface, which is explained by noting that 
our expression would fail as soon as the deformation of 
the contact lines became comparable to the distance be- 
tween them or to the drop size. In that case, the total 
capillary force might saturate, which supports the maxi- 
mal capillary force expression frequently used to describe 
detachment above a given drop size. 



CONCLUSION 

The wetted contact of a free interface with a cylindrical 
object was studied both numerically and experimentally. 
Quasistatic simulations give the deformation of the inter- 
face and the position of the contact line, in good agree- 
ment with experiments on a soap film in contact with a 
cylinder of liquid. Above a limiting velocity, the liquid jet 
applies a dynamical "total wetting" condition to the in- 
terface, with what we might call a dynamically-generated 
capillary force, which is also of paramount importance 



in explaining the refraction of the jet, as described in 
Kirstetter et al. [2l[ . From the contact line position, the 
force exerted by the interface on the object (and vice- 
versa) is calculated. We observe that whatever the in- 
clination of the cylinder, the force remains vertical and 
its amplitude is given by F^ = 2jttR/ cos{9i). This ex- 
pression should be useful in numerous problems in both 
fiuid and solid mechanics, and in future work it could be 
extended to study the partially wetted case, validated by 
the sort of simulations described here. 
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which gives the force components along the three axes of 
the reference frame of the cyhndcr TZ' {Ox'y'z'): 



F.'/l 



Fy'h 



sH0')%U0' 



-cos{e')^]de' 



F^,/j= j> Rde'. 



After some calculation, the contour integrals can be writ- 
ten in cylindrical or cartesian coordinates as: 

F^V7 = - <j> {z'cos{e'))de' 

z'd(sin(6i')) 



Fyh 



R 



-Rf/'y' 



{z' sin{e' ))d9' 



z'd{cosie')) 



z'dx' 



F^./j^rJ) d0' 



c 
2t:R. 



APPENDIX: EXPRESSION FOR THE FORCE 



Expanding the integral for the force (Eq. [S]) leads to 



F/7 



dz' 
Sr' A ( Ree' + -^e,> ] dO' 



dz' 



dO' 



Re, 



+ sinie')^e,,-cosie')^ey,ye\ 



The last expression for each component shows that F^i 
and Fy' are directly related to the projected areas of the 
contact line in the (Oy' z') and (0.t' z') planes respec- 
tively. With Ay'z' and Ax'z' the absolute values of these 
oriented areas, as defined previously, we have by symme- 
try (Fig. Is]) that Ax'z' = and /^ z'dy' < 0. This leads 
to the final expression for the force in the reference frame 
of the cylinder TZ' {Ox'y'z'): 

2j7rRfy''' 



F, 

Fy'^0 

Fz' = 277ri? 



27ri?2 



Writing the force components in the reference frame 
of the interface TZ can be done using the change-of- 
coordinates formulae (Eq. [I}, which then gives Eq. |6l 



